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DEFLECTION AND STRESS ANALYSIS OF THIN SOLID WINGS OF ARBITRARY PLAN FORM 
WITH PARTICULAR REFERENCE TO DELTA WINGS 3 


By MANUEL STEIN, J. EDWARD ANDERSON, and Jonn M. HEDGEPETH 


SUMMAR Y 


The structural analysis of arbitrary solid cantilever wings by 
small-deflection thin-plate theory is reduced to the solution of 
linear ordinary differential equations by the assumption that 
the chordwise deflections at any spanwise station may be ex- 
pressed in the form of a power series in which the coefficients are 
functions of the spanwise coordinate. If the series is limited to 
the first two and three terms (that is, of linear and parabolic 
chordwise deflections, respectively, are assumed), the differential 
equations for the coefficients are solved exactly for uniformly 
loaded solid delta wings of constant thickness and of symmetrical 
double-wedge airfoil section with constant thickness ratio. For 
cases for which exact solutions to the differential equations 
cannot be obtained, a numerical procedure is derived. Experi- 
mental deflection and stress data for constant-thickness delta- 
plate specumens of 45° and 60° sweep are presented and are 
found to compare favorably with the present theory. 


INTRODUCTION 


One of the present trends in the development of high-speed 
airplanes and missiles is toward the use of thin low-aspect- 
ratio wings. The structural analysis of these wings often 
cannot be based on beam theory since the structural defor- 
mations may vary considerably from those of a beam and, 
indeed, may more closely approach those of a plate. In 
cases where the wing construction is solid or nearly solid the 
use of plate theory in the analysis is particularly valid, and 
it is this type of wing which is considered in the present 
report. 

Exact solutions to the partial-differential equation of plate 
theory are not readily obtained, especially for plates of 
arbitrary shape and loading; however, a number of approxi- 
mate solutions to specific problems on cantilever plates have 
appeared in the literature (see, for example, refs. 1 to 7). 
Of the approaches used in these references, only the one in 
references 6 and 7 is readily applicable to plates of arbitrary 
plan form, thickness distribution, and load distribution; thus 
it is the most useful one for the analysis of actual wings. 

In reference 6 the cantilever-plate problem is simplified by 
the assumption that the deformations of the plate in the 
chordwise direction (parallel to the root) are linear. By 
minimizing the potential energy of the plate, the partial- 
differential equation of plate theory is replaced by two 


simultaneous ordinary differential equations for the spanwise 
variations of the bending deflection and twist. In reference 7 
the same ordinary differential equations are obtained in a 
different manner. Refinement of the analysis by inclusion 
of the effect of parabolic, cubic, or higher-order chordwise 
camber terms is indicated in reference 6, and as the order of 
refinement is increased a corresponding increase in the num- 
ber of ordinary differential equations is obtained. 

In the present report, which is an extension of reference 6, 
a general set of ordinary differential equations is presented 
which may be used to obtain any desired degree of approxi- 
mation to the deflection of the plate. These equations are 
solved exactly for several cases of delta plates under uniform 
load first by considering linear chordwise deformation only 
and second by including the effect of parabolic chordwise 
camber. Comparisons are drawn between the stresses and 
deflections computed from the equations of each approxi- 
mation and also with some experimental results. 

The differential equations presented contain coefficients 
that depend on the plan form and stiffness distribution of the 
plate and on the loading. In this report, the plates con- 
sidered in detail have coefficients such that the differential 
equations can be solved exactly; however, in cases for which 
exact solutions cannot be obtained a numerical procedure 
must be used. One such procedure is derived and its 
accuracy is demonstrated. 


SYMBOLS 

l length of plate measured perpendicular to root 
C root chord of plate 
p lateral load per unit area, positive in z- 

direction 
t local thickness of plate 
te average thickness of plate 

3 
D local flexural stiffness, itt 
D flexural stiffness based on average thickness, 
Ely 

12(1— p°) 
E modulus of elasticity of material 
u Poisson's ratio 
Ww deflection of plate, positive in z-direction 
wY coordinates defined in figure 1 


1 Supersedes NACA TN 2621, “Deflection and Stress Analysis of Thin Solid Wings of Arbitrary Plan Form With Particular Reference to Delta Wings” by Manuel Stein, J. Edward 


Anderson, and John M. Hedgepeth, 1952, 


2 REPORT 1131— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Kx, y) 


AOS 
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y=0 (xy _ 
bear feo 


FIGURE 1.—Coordinate system used in the present analysis for a canti- 
lever plate of arbitrary shape with arbitrary thickness variation. 


Pn function of x, coefficient in power series for 
. N 
deflection w= >} ¢,(x)y” 
n=0 


Cr(x), cala) functions defining plan form (see fig. 1) 


Ly variable obtained by transformation m=1—5 
Or, Oy normal stresses 
Tij shear stress 
o maximum principal stress 
l J2 
A aspect-ratio parameter, | 5 (1— u) 
RESULTS 


The derivation of the general set of ordinary differential 
equations is given in appendix A. The general procedure 
outlined in reference 6 is followed; that is, the deflection 
of the plate w is expanded into a power series in y the chord- 
wise coordinate with coefficients which are functions of æ 
the spanwise coordinate (see fig. 1) 


W= g(x) + or(z)y +e2l(x)y?+ ... +en(x)y® (1) 


Equation (1) is substituted into the expression for the po- 
tential energy of the plate-load combination which is in 
turn minimized by the calculus of variations with respect 
to each of the coefficients p,. The results of the variational 
procedure appear as N-+1 simultaneous differential equa- 
tions with the coefficients ¢, as unknowns. 

By taking a sufficient number of terms in the expansion 
of w, the resulting differential equations can be used to ob- 
tain any desired degree of accuracy in the solution for the 
deflections of any given cantilever plate subjected to an 
arbitrary lateral load. Of most interest, perhaps, are the 
particular cases for N=1 and N=2, which are obtained 
from the general set of equations and are simplified in 
appendix A. The case for N=1 (also derived in refs. 6 


and 7) includes linear chordwise deflections, and the case for 
N=2 takes into account parabolic chordwise curvature. 
Although for most practical problems the solution by the 
parabolic theory should be adequate, cases might exist in 
which cubic, quartic, or even higher-order chordwise terms 
should be included, depending on the convergence of the 
series for the configuration considered. 

The particular equations for N=1 and N=2 are used to 
determine the deflections and stresses of the following can- 
tilever plates subjected to uniform lateral load: 

(1) A 45° delta plate of uniform thickness 

(2) A 60° delta plate of uniform thickness 

(3) A 45° delta plate of symmetrical double-wedge airfoil 
section with constant thickness ratio 

Fortunately, for these configurations, the solution can be 
carried out exactly by both the linear and parabolic theories, 
and the details of these exact solutions are included in 
appendix B. In general, however, exact solutions cannot 
be obtained and some numerical method must be used. 
One such method, based on replacing derivatives by their 
first-order-approximation difference forms, is derived in 
appendix C. 

A summary of the results for the three particular problems 
is shown in figures 2 to 11. Deflections obtained by the 
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- Fraure 2.—Deflections of a 45° delta plate of uniform thickness 


uniform load. 
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linear theory and the parabolic theory for the three con- 
figurations are compared in figures 2, 3, and 4. Stresses 
obtained by the linear theory and the parabolic theory 
for the three configurations are compared in figures 5, 6, 
and 7. Where available, experimental deflections and 
stresses are also shown in these figures. The details of the 
procedure used to obtain the experimental deflections of 
the 45° and 60° uniform-thickness plates and the experi- 
mental stresses in the 45° uniform-thickness plate are con- 
tained in appendix D; whereas the experimental root stresses 
for the 60° uniform-thickness plate were obtained from 
reference 8. Figures 8 to 11 present the comparison be- 
tween deflections and stresses computed from the exact 
solutions of the differential equations and those computed 
from the numerical solutions of the same equations. 


DISCUSSION 


The results shown in figures 2 and 3 indicate that, with 
regard to deflections, either the linear theory or the parabolic 
theory is adequate for the case of a constant-thickness delta 
plate subjected to a uniform load, the comparison being some- 
what better for the 60° plate than for the 45° plate. If 
accurate slopes in the chordwise direction (angle of attack) 


Experiment | 


wD e 
} theory 


o—— — Linear 


i 5 | > Parabolic 
03 O i 
= 


Re) 


Figure 3.—Deflections of a 60° delta plate of uniform thickness under 
uniform load. 


are desired, however, the parabolic theory must be used 
because the error in the angle of attack as computed by the 
linear theory is as much as 30 percent (see figs. 2 and 3). 
The appreciable anticlastic curvature, evidenced by the 
experimental results of figures 2 and 3, may be important 
aerodynamically and is, of course, not taken into account by 
the linear theory. 

The apparent convergence of the aforementioned series 
in the case of the double-wedge-section plate (see fig. 4) 
implies that the linear theory is adequate for this case. The 
lack of chordwise curvature in the result obtained by the 
parabolic theory is attributable to the fact that the natural 
tendency of the plate to have anticlastic curvature is canceled 
by the opposite tendency of the thin edges to bend down 
under the load. Unfortunately, no experimental results 
are available for this configuration. 


In figure 4 the plate stiffness D in the nondimensional 


parameter wD/pl* is the local value of D at a point where 
the thickness is equal to the average thickness of the plate 
as a whole. Thus the results of figure 4 are comparable 
with the results of figure 2 on an equal-weight basis. It 
can be seen that the deflections of the double-wedge-section, 
constant-thickness-ratio plate are everywhere less than 


o-—- — Linear 


Porabolic }theory 


10 


We 


Figure 4.—Deflections of a 45° delta plate of symmetrical double- 
wedge airfoil section and constant thickness ratio under uniform load. 
D= Etas A 

12(1— p’) 
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those of the uniform-thickness plate although they increase 
rapidly near the tip. This curling-up or singularity in slope 
at the tip is a result of using a small-deflection theory and 
probably would not be so marked in an actual case. 

The stress results for the 45° and 60° uniform-thickness 
delta plates indicate that both the linear and the parabolic 
theories are adequate and that the parabolic theory is better 
than the linear theory only near the root. It should be noted 
that, although the maximum principal stress over a large 
part of the 45° plate is plotted in figure 5, only the stresses 
7=0.0087 of the 60° plate 
are plotted in figure 6 since only these stresses are given in 
reference 8. (The maximum principal stress and the stress 
normal to the root are theoretically equal at the root since 
the root shear stress is zero.) 

Experimental data are lacking for the double-wedge-section 
delta plate and, therefore, only theoretical stresses are shown 
in figure 7. As in the case of deflections, the results obtained 
from the linear theory and those obtained from the parabolic 


normal to the root along the line 


theory are almost coincident, the difference being greatest 


near the root. Figure 7 has also been plotted so that the 
results are directly comparable with those for the 45° 


20 


9 Experiment | 


—— — Linear 
Parabolic Theory 
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le 


FIGURE 5.—Maximum principal stresses in a 45° delta plate of uniform 
thickness under uniform load. 


uniform-thickness plate in figure 5 on an equal-weight 
basis. As can be expected, the double-wedge-section, 
constant-thickness-ratio plate is a better design structurally ; 
the stresses in the double-wedge-section plate are everywhere 
smaller and are almost constant in the spanwise direction. 

The theoretical results in figures 2 to 7 have been obtained 
from exact solutions of the differential equations of the 
linear and parabolic theories. In order to test the reliability 
of the numerical method derived in appendix C, the differ- 
ential equations were also solved numerically. The results 
shown in figures 8 and 9 indicate that the agreement is good 
between the numerical solution in which five equal intervals 
were used and the exact solution of the differential equations 
for the case of the 45° uniform-thickness plate. The same 
good agreement can be expected in other cases where the 
thickness and load distributions are not too erratic and where 
the plate stiffness does not go to zero at the tip—that is, 
when no singularities appear at the tip. 


2.4 


O Experiment (ref. 8) 
—— — Linear 
Parabolic 


} Theory 
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FIGURE 6.—Normal-stress distribution near the root (at z =0.0087 ) of 


a 60° delta plate of uniform thickness under uniform load. 
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Since the efficacy of the numerical method depends on 
how well parabolic arcs fit the various functions between 
stations, serious error can result from blind application. An 
example of the seriousness of these errors and of the manner 
in which they can be remedied is shown in figures 10 and 11. 
In these figures a comparison is made between exact and 
numerical results obtained on the 45° double-wedge-section, 
constant-thickness-ratio plate. As can be expected, the 
five-station numerical solution fails to follow the exact solu- 
tion in the neighborhood of the singularity at the tip. Since 
the region of trouble is localized at the tip, a reasonable 
remedy would be to decrease the spacing of the station points 
near the tip. This decrease in spacing may be accomplished 
either by using a greater number of equally spaced stations 
or by using unequally spaced stations crowded near the tip. 
The increase in accuracy obtained by increasing the number 
of equally spaced station points to ten is shown in figures 
10 and 11. 


CONCLUDING REMARKS 


The general method presented herein for finding deflec- 
tions and stresses of solid or nearly solid wings is, in principle, 


2.0 


—— — Linear } 
Parabolic Theory 


NU 
ann 


FIGURE 7.—Maximum principal ee in a 45° delta plate of symmet- 
rical double-wedge airfoil section and constant thickness ratio under 
uniform load. 


capable of yielding arbitrarily accurate results for any con- 
figuration. It is seen that, for the examples considered, 
only the first two or three terms in the series expansion need 
be considered to obtain adequate accuracy. In addition, 
for most practical plate-like wings with clamped roots the 
first two or three terms will probably be adequate, although 
problems may exist wherein more terms are needed. 

The numerical procedure, derived for application in cases 
where exact solutions cannot be obtained, gives good agree- 
ment when compared with exact solutions if enough stations 
are taken along the span. The necessary number of stations 
is dependent on the type of thickness and loading distribution 
considered, five equally spaced stations being enough for 
the uniform-thickness delta wing subjected to uniform 
loading and ten being necessary for the double-wedge-section, 
constant-thickness-ratio delta wing subjected to uniform 
loading. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LANGLEY FreLD, Va., November 30, 1951. 
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Figure 8.—-Numerical and exact solutions of the differential equations 
for the deflections of the free edges of a 45° delta plate of uniform 
thickness under uniform load. 
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N 
(obtained by assuming linear chordwise deflections) for the maximum 0 ý Y 
principal stresses along the free edges of a 45° delta plate of uniform 
thickness under uniform load. 
“Ol 2 4 


a a 

Figure 9.—Numerical and exact solutions of the differential equations 
6 8 I.O 

x/l 

FIGURE 10.—Numerical and exact solutions of the differential equations 
(obtained by assuming linear chordwise deflections) for the deflee- 
tions along the free edges of a 45° delta plate of symmetrical double- 
wedge airfoil section and constant thickness ratio under uniform load 
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FIGURE 11.—Numerical and exact solutions of the differential equa- 


tions (obtained by assuming linear chordwise deflections) for the 
maximum principal stress along the line IQ) of a 45° delta 


plate of symmetrical double-wedge airfoil section and constant 
thickness ratio under uniform load. 


APPENDIX A 


DERIVATION OF DIFFERENTIAL EQUATIONS © 


The structure considered herein is a thin, elastic, isotropic, 
cantilever plate of arbitrary plan form and slowly varying 
thickness subjected to distributed lateral load (see fig. 1). 
By assuming that the deflection of the plate can be repre- 
sented by a power series in the chordwise coordinate and 
by applying the minimum-potential-energy principle, a set 
of ordinary differential equations in the spanwise coordinate 
is obtained from which the coefficients of the power series 
may be determined. First the general set of equations is 
derived; then the particular equations for the cases of linear 
chordwise deflections and parabolic chordwise deflections are 
deduced from the general set and simplified. 


266368—-54——-2 


General equations. —The potential energy of the system 


under consideration is 


| a? Co (£) 2 a 
Potential energy = | J ) ; k id (E) + 2) 6 (Sr) + 
l 0 y CI 


a paro) aww dy da (AD 


in which A 
ix, yl? 
DyN 


and p(x,y) is the distributed lateral load. 
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The assumption is made that the deflection w can be represented by the power scries 
N 
w= pa Gn(x)y” (A2) 
n= 


Substitution of this expression for w into equation (A1) gives 


`Z 1 N N 
Potential energy= | deba 2 2 [Amyn pm Pr Hmn (m — 1) (1) AUm+n—3 Cm PaT 
0 m=0 n= 
N 
2u4n (n— anini m Ca +2(1 ES PIMNAm+n—18m En | z2 Pasien? (A3) 
n= 


in which 
Ca (x) 
a.—| : D(x yyy dy C= 5-24 2N +1) 


cy (r) 


Ca (1) (A4) 
p= fo p(x, yy’ dy (r=1,2,...N+1) 


and the primes denote differentiation with respect to x. 
Minimization of the potential energy by means of the calculus of variations gives 


¿(Potential energy) =0 


L 


l 1 N N 
-=f dx 5 27 2o [Omang (Om Spn Fenr pm + mn (M—1) (N= 1) Aine n—3(PmdGn F Pnp m + 


m=0 n=0 


N 
2un(n— Dam+n-1(Cm 06m IÓ T 2(1 ES p)MNAm+n— Cm 0Pn + Gn’ bOm’)| = >i Putt 5) 


n=0 


Integrating by parts and collecting terms results in 


l N N 
= f, á 2 On = (Am+n+10m ) +m (m—1) (amtn-1Pm) = 2(1—u)mn (Amtn-18m) Hun (n— 1) Amtn-19m F 
n= m= 


= 
y 


N ¿ N N 
mn (m — 1) (n— 1) Um+n—3 enl—Pnss} + p> DPn PR [Amyn Pm tum (m — 1) Amin) — > OPn oa A n ITF 
n= n= m= 


m=0 
| l 
um (m 4 1) CA Gy) = 2 (1 = p) MNAmn+n—1 n't (A5) 
0) 
Everywhere in the region of the plate, except at the boundary x=0, the variation of w is arbitrary. At x=0 the 
cantilever boundary conditions 


Ow 


=>; 70 


W 


yield 
Pn(0)= px (0)=0 (n=0, Lice tie N) (A6) 


and therefore the variation in these quantities must also be zero at x=0. 
Equation (A5) is then satisfied if, in addition to equation (A6), 


N l 
27 [(Am+n+10m )” + um(m— 1) (Unida) 2(1 oe p)MN(A m+n—18m.) + un(n— Dainim T 


m=0 
mn(m—1)M— Da n+r-30n1=Pr+1 (n=0,1,...N) (A7) 
D lansarien" HMM — Da n+a-1 mbr =0 (m=0,1,...N) | (A8) 
and 
D [ansar en Y + umn 1) (Am+n—1€m) —2(1 — MBA m1 lz-1=0 (n=0,1,...N) (A9) 


Equations (A7) form a set of N+1 simultaneous ordinary differential equations for the functions ¢,(x). The functions g, 
are completely determined by these differential equations and the boundary conditions (A6), (A8), and (A9). 
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Particular case of linear chordwise deflections.—If 
N=1, the deflection function becomes 


W=P + Y (A10) 


a linear function in the chordwise direction, where g is the 
bending deflection and e, is the twist. Equations (A7) 
become 


(ap + (asp) =p: 


(apo) + (ap) —2(1— u) (0191) =p 


(A11) 
(A12) 


The root boundary conditions, given by equation (A6), 
become | 

90(0) = go’ (0) =p (0) =p," (0) =0 (A13) 
The tip boundary conditions, given by equations (A8) and 
(A9), become 


(ap "+21? 7-10 (A14) 
(dogo +. d3¢1/") 1 (A15) 
Largo") + (mg) ]z==0 (A16) 
[(d2qo"")’ +(a30 Y —2(1—wa pl, .=0 (A17) 


Equations (A11) to (A17) are the differential equations and 
corresponding boundary conditions presented in reference 
6 Gf only distributed load is considered) where the symbols 
W and 0 are used instead of g and gı, respectively. 

If equation (A11) is integrated twice and the boundary 
conditions (A14) and (A16) are used, 


if Qe ff 1 de d 2 
A 
Po ay pı a), oa 


Substitution of pọ” into equations (A12), (A15), and (A17) 
gives 


l i ff 
(bier! 20.0) mp =p (2 f f j da) (A19) 


(A18) 


(bi¢1"")--1=0 (A20) 
[(01p1)—2(1— war ls. .=0 (A21) 
in which 
2 
b, 4 


If equation (A19) is integrated once and the boundary 
condition (A21) is used, 


l l l A 
Boy 20) == | Pa dE fip da*) (A22) 


The differential equation (A22) is a second-order differential 
equation in y”. The twist e, and then the bending deflection 
po are obtained by solving equations (A22) and (A18), respec- 
tively, by applying the boundary conditions (A13) and (A20). 

Particular case of parabolic chordwise deflections.—The 
effect of parabolic chordwise camber may be included by 


letting N=2 in the general power series (eq. (A2)). If N=2, 


the deflection function becomes 
w= po H Yp +Y*02 


Here pz represents the spanwise distribution of parabolic 
chordwise camber. For this case the differential equations 
(A7) become 


(Aigo )” + (ap) + (d3p2”")” +2u(41p2)” —P1 (A23) 
(a2p0")"+ (ap) + (dap2”)” 2u (47p2)” — 
2(1 — p)[(di¢1')’ +2(a2p2)”] =p, (A24) 


(A3po”) vt. (ap) + (Aspo”) "+2ula1p0” + Arp,” + azp” + 
(A3p2)"] 4 (1— pw) [(degy’)’ +2 (d3p”) ‘1+ 4a g == P3 
with the boundary conditions 


po(0) =p (0) = ¢1(0) = 91’ (0) =p.(0) =p," (0)=0 (A26) 


(A25) 


(Aigo" + d2g1” + A3¢2” + 2140182) 2:=0 (A27) 
(azp0” asp” Hup” + 2ud2¢2) 2-1-=0 (A28) 
(A3p0” + asp” +4502” + 2yd3¢2)2-1=0 (A29) 
[(d1po”) + (dagr”)’ + (asp) + 2u(Gigr)]2-1==0 (A30) 

(apo) + (ap) + (tp: y +2u(47p,)' — 
2(1—u) digi’ +242p2)l7-:=0 (A31) 

[(azp") + (ap) + (asp) +2u(a3p.)' — 
4(1—u) (ap: + 24362 )+-.=0 (A32) 


If equation (A23) is integrated twice and the boundary 
conditions (A27) and (A30) are used, 


T 2 Y E ” I ae hs 2 ‘ 

Po =—— oO —— gy —2up tH — pı dx (A33) 
ay ay) di Jue Jr 

Substitution of py” into the remaining differential equations 


and boundary conditions results in 


(010 + (b2g2")” — 2(1 — u) Kargi y + 2(a2@’)’] 


a (it? n 
= Pa— aff Pı dæ?) (A34) 
Qi Jr dx 


(bapi Y” HO Y” — 4(1 — u) [(arp Y + 2(a3@0’)’] +41 — pare. 


one [ifin aee f Ln] ase 

(bipi bp” )r-1=0 (A36) 
[ry + Op") 20a) mg 420002 l:=0 (A37) 
(dapi” Hbg” )r=1=0 (A38) 

eriten aa a A 


gı (0) =¢1' (0)=p(0)= p (0) =0 (A40) 
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in which 


d 

b= E 
ay 

oh 

bo = d4 ee 

ay 

2 
a3 

b, =a; — — 
3 5 a 


If equation (A34) is integrated and the boundary condition 
(A37) is used, 


(dipi”)' +b” Y — 2(1 =h) (aigi + 24292’) 


; l a Ll 7 
== Par (E f f mde) (aan 
T 1 z x 
Thus q, and e, are obtained by solving equations (A35) and 
(A41) with the boundary conditions (A86), (A38), (A39), and 


(A40). Subsequently, g can be obtained by solving equa- 
tion (A33) with the boundary conditions g(0)=¢’ (0) =0. 


Stresses.—After the approximate deflection of the plate 
is determined from equations (A18) and (A22) or from equa- 
tions (A33), (A35), and (A41), the extreme-fiber stresses may 
be calculated from the well-known equations of thin-plate 
theory, which are (see, for example, ref. 9): 


_6D (Vw, Mw 
ae Nee Oy 


6D (òw DW 
Oy E dy? Pay? 


BAD dw 
A e Ox dy 


The maximum principal stress ¢ at any point in the plate 
can be determined from 


: l ~— 
pees ae i 9 V(G2— oy +4 72 


APPENDIX B 


EXACT SOLUTIONS OF DIFFERENTIAL EQUATIONS FOR SOME SPECIFIC DELTA-PLATE PROBLEMS 


The differential equations of appendix A for linear and 
parabolic chordwise deflections are solved exactly for uni- 
formly loaded delta plates of constant thickness and of sym- 
metrical double-wedge airfoil section with constant thickness 
ratio. The equations for deflections obtained by the linear 
theory are presented in terms of the aspect-ratio parameter 
A for both kinds of delta plates. The equations for deflec- 
tions obtained by the parabolic theory are presented for 


y v3 


a. and “a7 with => for the constant-thickness delta plate 


; 1 
and for : =1, also with u= for the delta plate of symmet- 


rical double-wedge airfoil section with constant thickness 
ratio. 

If the z-axis is passed through the edge perpendicular to 
the root and the substitution = 1 =5 is made, the differen- 
tial equations are clearly of the homogeneous type for which 
the solutions are of the form 2,7, where y is a constant. For 
the configurations considered, the functions that define the 
plan form (see fig. 1) are then c (x)=0 and cz(x) =cx1,, where e 
is the root chord. In all the equations of this appendix the 
primes denote differentiation with respect to the new inde- 
pendent variable 2. 


DELTA PLATE OF UNIFORM THICKNESS UNDER UNIFORM LOAD 


Since the stiffness D is a constant for uniform-thickness 
plates, the coefficients in the differential equations (sce eq. 


(A4)) become 


¿DE a dl 
by =a; a 2 (B1b) 
b3= ite LP (Bid) 

pre ay (Ble) 


Solution for linear chordwise defiections.—If the co- 
efficients given by equations (B1) are substituted into equa- 
tions (A22) and (A18) and the independent variable is 
dí 


l 


changed to zı=1—7» the following equations for linear 


chordwise deflections result: 


(xg Y — 16° zig = — 2 Dc Pe (B2) 
4 
go = -5 ne” + Pr Di (B3) 


where 7 
a ae 
ar NE (1 —u) 


The boundary conditions to be used with these equations are 
obtained from equations (A13) and (A20) and are 


gll = g MD=a)=e(1)=0 (B4) 
(oa 0 B5) 
The general solution of equation (B2) is 
/ = are ae l 
gy’ = Ayr? !4+ Aor, 77 AS a (B6) 


where 
y= V1 16% 


and 4, and 4 are arbitrary constants. Since A? is inherently 
positive, the boundary condition (B5) requires that A,=0. 
One integration of equation (B6) and the application of the 
conditions ¢;(1)=¢’(1) =0 yields 


7 1 plt fayy—-1 2-1 
e l 
If equation (B3) is solved for g) with the conditions 
go(1) = g (1) =0, the result is 
a pli 1l 2 rR Ay? ( at > 
0" 8D 1-2 É ON OE E” 
yY 1 1 i) 
“— ( 1—97,— B8 
Y EE ae Gi 


Substitution of equations (B7) and (B8) into the equation 
w= Got Yo 


gives the expression for the deflection w of the plate under 
the assumption of linear chordwise deflections. 

Solution for parabolic chordwise deflections.—If the co- 
efficients given by equations (B1) are substituted into equa- 
tions (A41), (A35), and (A33) and the independent variable 
is again changed to 1 =1-7 the following equations for 

11 
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parabolic chordwise deflections result: 
INT e r E [t € 
(10! ) +21 tp?! — 1 6? (zipi +xtco, ) — 2 e e 
(B9) 


tt 16 Z a : / 4 E WN! 
(xip) tis (21°C gy” JP 16M (tor) +a (x "ego ) |+ 


1 2unr? 
ae i wet a e De z? (B10) 
ff C 2 pl* 
E Lipi” - rep? ~ Qul CS y (B11) 


The boundary conditions to be used with these equations are 
go(1)= go =A) (1)=e(1)=p'(1)=0 (B12) 


(113 91" + 21 4e¢2"")7,-0= 0 (B13) 


tf 16 21 
(ata +È ni'eo: ) _=0 (B14) 


ff 4 
[Cos pr) + (a cos) — 16y? (a? Pı ole step, ) Fa 
(B15) 


The homogeneous solutions of the simultaneous equations 
(B9) and (B10) are of the form 


gy =Ax Y? 


go= Bay 


Substitution of these expressions into the homogeneous parts — 


of equations (B9) and (B10) leads to the following charac- 
teristic equation from which \ may be determined: 


76+ 16N)y4| 320 (44,54) M+48007+9 [y 


4 | 1280 57 ce r+ 80 (44,74) MH 96X24 1 |o 
(B16) 


and gives the following relationship between A and B: 


Fie TA y AAA aL Jez 


The particular solutions for uniform loading are given by 


o =A,2? 
p.=bB,x1* 
where 
4 SELI uio 2 ce ; 
Jozi = “1 pi 
Á. pr : 
g E ont Di mane i 
n 2 C~- SS f 
B,=—}- pl 


TEN "oy 1) (BN aN nee 


The general solution is the sum of the homogencous solu- 
tions and the particular integral 


6 
Pı = Apt 14+ Api 
n= 
6 
p => Bt I+ Bz? (B17) 
n=1 


where the values y, are the six roots of the characteristic 
equation (B16) and the coefficients A, and Bb, are the co- 


efficients corresponding to each of these roots. After inte- 
gration gı becomes 
9 Lyn ry A 
n= n 


The general solution for gy from equation (B11) is found 
to be 


6 
p= Cam PL Ozi C21 + O, (B 19) 
where, for n=1, 2,. .. 6, 
SEE Ey k 9 a An+$| m6 2+7 = la, ! 


and 


= 20 2uN pi 
Les 5| 4+ (+E ~) By A 


The coefficients A; to As, A, C, and C, must be determined 
by the boundary conditione (B12) to (B15). 
A complete set of coefficients is given in the following table 


for delta plates with Poisson’s ratio y equal to 1 /3 and with 
peasy and 42. Deflection curves plotted from these 


results are shown in figures 2 and 3 in which the 45° plate 


l l 3. 
corresponds to —=1 and the 60° plate corresponds to o 
C 
De De D 
Ym Aom pú Bm pl Con 
m fe A e —-— =e- e-m ene t ee name il ae mena oe o A ——— Km — e eae aaa 
ker NS A=1 pr 13 a 1 V3 A=] p= V3 
3 3 3 3 
1 | 2.7034 | 1.5671 | 0.7378 | 0.09632 |—0.3133 |—0.1022 |—0.02931 |—0.003223 | 
2 4. 9437 3. 6347 . 02411 . 3707 —, 03039 | —. 4313 . 003074 . 01347 j 
3 | 8,3816 | 4.7258 | .03827 |—.1766 | —.006298| 107379 | 1000486 | | 909817 
4 |—2. 7034 |—1. 5671 | 0 0 0 0 0 0 
5 |—4. 9437 |—3.6347 | 0 0 0 0 0 0 
6 |—8. 3816 —4. 7258 | 0 0 0 0 0 0 
rN ai eee —.8000 |—. 2903 -3500 . 4597 04167 -004032 | 
70 OO NINAS —.01557 |—.92924 | _..... | —.07152 | ~108354 | 
PG case ae AECA WW ealiismins la TE heen . 05668 E 06692 
r Mas MERA a 


Substitution of equations (B17), (B18), and (B19) into the 
equation 
W=gotyaty es 


gives the expression for the deflection w of the plate under the 
assumption of parabolic chordwise deflection. 
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DELTA PLATE OF SYMMETRICAL DOUBLE-WEDGE AIRFOIL SECTION WITH 
CONSTANT THICKNESS RATIO UNDER UNIFORM LOAD 


For a delta plate of symmetrical double-wedee airfoil 
section with a constant thickness ratio the thickness is a 
function of z and y and is given by the following equations: 


(0 <y< 3) 
q CX 
j=. (2—2) (S sysen) 


where ta is the average thickness. From these expressions 
for the thickness the stiffness can be found and the coeffi- 


i= Bla 2 
C 


(A18) for lincar chordwise deflections may be solved for gı 
and po. The steps in the solution are the same in form as 
those for the uniform-thickness plate and the resulting 
equations are 


(B21) 


cients in the differential equations become r 
ne ys 
La 4 A] a 
SER 2 (* 
27 De? A i 973 (21 log, 21 + 1 — 21) pr (B22) 
g= > Fi ne D 
8 a 
9De , where 
y == ry PE ~ 
5 5y 2 
B v=/(3) + 80° 
siDet ., 
Eon ¡ ¡ 
Solution for parabolic chordwise deflections.—By use of 
2673 De? the coefficients given by equations (B20), equations (A41) 
o (B20) | | 
5= 4480 (A35), and (A33) for parabolic chordwise deflections may be 
a solved for gı, ¢2, and gp. The steps in the solution are 
E pi again the same in form as those for the uniform-thickness 
80 plate and the resulting general expressions for gı, 2, and 
A , Po- are 
a 380 “i 6 r Yn — l 
— pensa m2 É 
y 2613D a a= 25 A, +4, log. z, (B23) 
¿7 22400“! E 
C” 
p,=* > >. 1 
N => Dar” 2+6, z (B24) 
n=1 1 
Solution for linear chordwise deflections.—By use of the a 
coefficients given by equations (B20), equations (A22) and a=} Crei" 2+0,2%; log. Oi +0, (B25) 
where the exponents y, are the roots of the characteristic equation 
, 25 871/ , 25 2 T a 320 + 01 2.225 a 
(rs) ada 2 Nr or (rom) + 
7/, 25 2 
For n=1, 2,...6, An, Bn, and C, are related by 
y OS 


B,= 


5 on 


Yar 


| y 
ND 


vale 


G LO DG) sor]? 
Antes ¢ | (7, 3) (1 Ji P 
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For uniform load the coefficients in the particular integrals of equations (B23), (B24), and (B25) are 


10 pit =F +10240°4+320 E M16 (547, x) OH) 


~ 27 De (2002-41) Croa T EA Mt) — 120 (162-41) (1002-41) 


2 2 Zur 
40 plt 15(16141)—2(202+ 1) (5+ >) 


27 De 2041) ($4 4 1024124320 a 120 (16024 1) (LON 1) 


p 


2 pl” 


O) =2 A E ez ee By toy ID 


The coefficients A, to As, A,, C,, and C, are again determined by the boundary conditions (A26), (A36), (A38), and (A39) in 
which the coefficients given by equations (B20) are substituted. 


For Poisson’s ratio u equal to 5 and pain], the solution of the characteristic equation (B26) leads to two real values and 
two pairs of complex conjugate values for y. The identity | 
r= =x cos (b log, x1) +12," sin (b log, x,) 
was therefore used to transform the terms involving the complex conjugate values into real form. If a and => the 


solution is 


Po = pt 0.0040702,° °°" — 0.00436325°" cos (2.825 log, xı) +0.0068932,2° sin (2.825 log, 21) + 0.000294 A 
e? 1 
4 

E O DO 0,0021342 07 cos (2.825 log, 11) —0.0063817,2:75 sin (2.825 log, x1) +- 0.01794 log, 21+ 0.001768 | 
C 


p= -5 |o. 00077 152°” —0.00007082,"°"* cos (2.825 log, 21) +0.0012342, sin (2.825 log, 21) + 


0.033312, log, zı —0.04096:x, + 0.04026 | 


APPENDIX C 
NUMERICAL PROCEDURE FOR SOLVING DIFFERENTIAL EQUATIONS 


In cases where the equations of the present theory cannot 
be solved exactly a numerical method must be used. In 
this appendix, equations (A19) and equations (A34) and 
(A35) are set up in difference form for numerical solution. 
Initially the assumption is made that the functions involved 
in the differential equations are continuous and nonsingular. 
In this case, first and second derivatives can be expressed by 
the standard difference forms 


CYN Yni 2Yn t+ Yn-t 
Ie 


Ee 


dy Eo Yn-; 
ea 


€ 


where e is the distance between equally spaced station points. 

In the following development five equally spaced span- 
wise stations are used; however, the extension to a different 
number of stations can be readily made. 


The quantity (big, y at the half-stations can be expressed 


po | 2 1 
pr”; 1 —2 
p” =; 1 
r 
O 3 
yf 
Yr 4 


First, consider equation (A19) resulting from the linear 
theory 


1 l US 
My 210 (ES f po dr =a 


Because of the nature of the tip boundary conditions for 
this equation, 1t can be conveniently put in the form 


where 
T= (b,¢,")’ —-201 —u)a pr” 


In finding the difference equation equivalent to equation 
(C1), the quantity (61017) is found in matrix form; from 
this expression is subtracted the matrix equivalent of 
2(1—y)di¢i’; the resulting expression for T is multiplied by 
a differentiating matrix; and the product is equated to the 
right-hand side. 


in matrix form as follows: 


P1-1 
1 P10 
md 1 Pır (C2) 
| o g 1 P12 


Pla 
€15 


where the second subscript denotes the station point, the ADO at the root station being 0 and at the tip 5. The 


root boundary conditions are now applied; namely, 


pı (0) = U= pio 
, —Q0—fu 1-1 
pı (0)=0 De 
Thus, after the values of gj.=0 and gı-1= gn are substituted, equation (C2) becomes 

eo 2 P11 
p1” mA | 1 12 
pr” e L *=2 1 $13 
er” | = 1 P14 
pra 1 Y 1 P15 
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Therefore, 
(0101 ”)o bio 2 Cn 
biyi”) by E 1 Pr 
(O1¢1")2 => by b =2 1 P13 (C3) 
(b1¢1")s bis LZ 1 P14 
(bigi Ja by, 1-2 1 ee 
One of the tip boundary conditions is 
(b101" r= ¿=0= (0,617); 
Thus, : 
(big, aye aa. I (brp1”)o 
(b1pp1”) 3,2 al 1 (b1p,”), 
(big1”) 5/2 =i —1 1 | (di¢1”)2 (C4) 
(dp) 7/2 sl 1 || ig”) 
(drpr”)oy2 — 1 | biy” Ja 
The matrix equivalent for the second term of T is 
(Qigy’)ip 41,17 1 P11 
(d191") 372 Q1,3/2 =l 1 P12 
2(1—H) (agi) sp UA ai, 5/2 | a. l $13 
(igi rz 4, 7/2 me l Pta 
(igi 72 1, 9/2 aul 1J | ¢15 
Therefore T becomes 
Tiz — 1 1 bio 2 
Ts 1 oh l bi i 1 
Tol =(= —1 1 biz e 3] — 
Tin == 1 bis 1 =2 l 
Topo — I bi, 1—2 1 
li, 1/2 l Pr 
2(1—p) Qy, 3/2 1 l P12 
oa = hy 5/2 —1 1 P13 (C5) 
47 s=] 1 Pla 
U, 9/2 sl 1_}) |215 


The right-hand side of equation (C1) can now be equated to the derivative of equation (C5); thus 


? 


du —1 1 Th) 
Yi e al 1 Ts 
413 e =A Ts 
i4 1 1J Trg 

Top 

In order to obtain qıs, the boundary condition 
T=0 
at x=/ must be used. In other words, T goes from T,» at station E to 0 at station 5. A straight line drawn between 


2 
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these two points would have the slope e. The value of qı; is considered to be this slope; therefore, 
Gu de: 1 Tig 
diz =l 1 Tp 
413 =" ak l Tp 
ds = 1 To 
(5/2 —11 Tog 
or 
qi 1 —2 1 Dro 2 
412 1 —2 1 bii —2 1 
M3 |= > 1-2 1 biz 1 —2 1 — 
da 1 —2 bis 1-2 1 
15/2 1 bi 1-2 1 
—1 1 11/2 l en 
—1 1 3/2 mal 1 P12 
2(— p) 
D = d Ay. 5/2 a a | P13 
— 1 1 dz — 1 1 Pra 
ee | Ai, 9/2 d 1 P15 


If the matrix multiplication is carried out, the difference equivalent of equation (C1) finally becomes 


qu Pr 
412 P12 
dis = e ¡0740 iD} P13 (C6) 
Gu Pta 
Q15/2 P15 
where 
2619+ 4b +b —2bu— 2b: biz 
— 2b — 2b  bu+t4brtbg  — 2b: — 2b biz 
[C]= bo —2b12—2b13  dio+4big+bi4S —2b13—2bj4 bu 
biz —2bi3— 2014 biz +4614 — 2b; 
Dr, — 2044 Dis 
— 1,12 H, 3/2 dy, 3/2 
Ur, 3/2 — a1, 3/2 l1, 5/2 1, 5/2 
[DiJ= A1, 572 — 1,51 U1, 772 1, 7/2 
01,7/2 — 4 1,7/2— QU, 972 M1, 9/2 


1,9/2 — A1, 9/2 


In order to determine g from gı, use must be made of equation (A18) 


tr 2 2 ad 
Po idz (1 
Qı z z p aA 


or, by use of the boundary condition (0) =’ (0) =0, 


x x 1 l l ip x (l: 
a= | f ~{ f | i} - gi’ da? (C7) 
0J0 Aids Jz ovo Ay . 
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In matrix form equation (C7) becomes 


Por 1 1/2 1/41 1 1 1 1 2771 1 1 1 1Tlo, 
Poz 1 1 1/2 1 1/41; 1 1 1 1 1 1 1 1l|p» 
gosi=et| 1 1 1 V2). 4 1/ay2 1 1 1 i i thingy |= 
Pos 1111 1/2 1 1 1 1/a13 1 1 1 1 lp, 
Los 1 1 1 1 1J[11/2 1 1 1 1 1/41, 1 1 p/2 
ÂE! 1/2 (a/a) 2 Or 
1 1 1/2 1 (a/a,); —2 1 Pro 
E: a 1/2 1 1 (42/41)» 1-2 1 Pr (C8) 
1 1 1 1 19 1 11 (42/43 Leo: j Ora 
1 1 1 1 1](1/2 1 1 1 1 (a/ai); 1 —2  1]les 


Thus, if the values of q, (which can be determined numerically or analytically according to preference and feasibility) are 
known, the values of p, can be found by solving equation (C6) and the values of gp in turn by means of equation (C8). 

The foregoing development applies to the case where only linear chordwise deformations are allowed. A similar pro- 
cedure is followed in expressing the differential equations pertaining to the parabolic theory in difference form; only the 
results are shown herein. 


The matrix equivalent to equations (A384) and (A35) is 


du | Pr 
(12 | Cr 
ie | C= RE ies 
Y ! P14 
1 A: Aa ae (C9) 
q21 $21 
22 | P22 
as | [att wa | Fie-8S >” [lon 
das (13) E] Pos 
q25/2 | | | P25 
where 
209 +4 Bit On2 = —2bn1 —2dnz bnz 
—20m1—20n2 Ont + 4bno-+Hbns —2bn2—2br3 bn3 
C= bn —2bm2—2bn3 bno tA Dag tbn —2br3—2bw4 bn 
bn3 —2bm3 — 20 m4 bngt+4bn4 —2bn4 
Bn —2bns Dis 
— En, 172 — Un, 3/2 Un, 3/2 
An, 3/2 — Un, 3/2 ln, 512 Un, 5/2 
[D,]= On, 5/2 An, 5/2 Un, 7/2 On, 7/2 
On, 7/2 — Gn, 772— En, 972 An, op 


On, 9/2 — An, 9/2 
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411 
Cra 
[1] = 413 
hig 
Qis 
and qı and q: are the right-hand sides of equations (A34) and (A35), respectively; that is, 
am (2 fi fine) 
Ll ae (ET ” 
pome finas (2 fine) 
zJ li Ja Jz 
With y, and g: known, py can be obtained by use of equation (A33) 
ee) id 1/2 1/aro 1 1 1 1 TIL 1 1 1 pu 
Po? 1 1 2 1 L/ay 1 1 1 1 1 1 1 Lf lpr 
w ltda 1/2 1 1 é 1 /ays 1 1 1 1 1 dilo. le 
ou. hac Se a 1/2 1 1 1 1/ay3 11 1 1] lp. 
ge A ALA a- i 1/ays 1 1] |p,5/2 
(d2/1)o 2 Pil 
(47/01): —2 P12 
(d2/d1)2 L =2 1 $13] — 
(a7/41)3 1-2 1 P14 
(2/a)4 1 —2 11 [015 
(a3/a1)0 2 $21 21 
(a3/a1), 2 1 P22 P22 
(d3/a1)2 | a 1 Pa| — 2 ue” | paz (C10) 
(a3/d1)3 LZ 1 P24 P24 
L (93/41), Le 11 [0 P25 


It should be noted that, as can be expected, the matrix 
equations (C6) and (C8) are merely special cases of equations 
(C9) and (C10), respectively. In addition, the square 
matrices in equations (C6) and (C9) are symmetric, a result 
that is consistent with the fact that the differential equations 
under consideration are self-adjoint. 

In the beginning of this appendix the assumption was made 
that the functions involved in the differential equations are 
continuous and nonsingular. The difference solution, how- 
ever, may be adequate for some cases in which this assump- 
tion is not strictly correct. For instance, the deflections of 
a plate with a discontinuous stiffness distribution could con- 
ceivably be not very different from the deflections of a plate 


with a continuous stiffness distribution closely approximating 
the discontinuous distribution except in the neighborhood of 
the discontinuity. The results yielded by the difference solu- 
tion in this case would be those associated with the continuous 
stiffness distribution. The number of stations may have to 
be increased, however, in order to minimize the inaccuracy 
introduced by the discontinuity or, in other cases, by a 
singularity. The case of the symmetrical double-wedge air- 
foil section, constant-thickness-ratio delta plate, discussed 
in the body of this report, is an example of a treatment of a 
singularity. In this case, although the solution is singular, 
adequate accuracy is obtained by the difference solution if 
ten equal intervals are used. 


APPENDIX D 
DEFLECTION AND STRESS EXPERIMENTS ON SOME TRIANGULAR CANTILEVER PLATES 


Test specimens.—The specimens tested were: (1) a 45% 
right-triangular plate clamped along one leg and (2) a 602 
right-triangular plate clamped along the longer leg. Each 
specimen, cut from 24S—T4 aluminum-alloy sheet of 0.250- 
inch thickness, had a length perpendicular to the clamped 
edge of 30 inches. 

Method of testing.—Figure 12, a photograph of the test 
setup, shows the methods of clamping, loading, and measure- 
ment of deflections. A 1,000,000-pound clamping load (held 
constant during the test) was applied to the root area of each 
specimen and a uniform load of 0.204 psi was applied by 
2-inch washers giving a tip deflection in each case of approx- 
imately 2 inch. | 

The deflections were measured by dial gages placed at the 
points indicated in figures 2 and 3. 

‘Stresses were obtained from the 45° specimen only. On 
this specimen, 13 resistance-wire rosette strain gages were 
placed at the points indicated in figure 5. The plate was 
loaded with 2-inch washers in four increments of 0.0847 psi 
per increment and the maximum tip deflection was 1.13 
inches. Readings of all the strain gages were recorded at 
each increment of loading. 

Analysis and discussion of data.—The deflection w was 
plotted in figures 2 and 3 in terms of the nondimensional 


Figure 12.— Deflection test setup of the 45° delta plate under uniform 
load. 
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parameter wD/pl*, in which the elastic constants were taken 


as £=10.6X10% psi and wy It was found that the dial- 


gage forces reduced the tip deflection of the plate by approx- 
imately 2 percent; however, since this error is of the same 
order of magnitude as that in the material properties and 
from other sources, no corrections are made in the results 
presented. | 

The readings of each of the 39 individual strain gages were 
plotted against load, and the slope of each of the resulting 
linear curves was taken as the average strain per unit load 
of the individual gage. The principal stresses were then 
calculated and plotted in figure 5 in terms of the nondimen- 
sional parameter ot?/pl’. 
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